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EXISTENCE AND STABILITY OF A SCREW DISLOCATION UNDER 

ANTI-PLANE DEFORMATION 

T. HUDSON AND C. ORTNER 

Abstract. We formulate a variational model for a geometrically necessary screw dislocation 
in an anti-plane lattice model at zero temperature. Invariance of the energy functional under 
lattice symmetries renders the problem non-coercive. Nevertheless, by establishing coercivity 
with respect to the elastic strain and a concentration compactness principle, we prove existence 
of a global energy minimiser and thus demonstrate that dislocations are globally stable equilibria 
within our model. 



1. Introduction 



m 

o 

^ Dislocations are line defects in crystalline solids which can be described by discontinuous 

^ displacements of a homogeneous crystal: to obtain the simplest forms of dislocations, the 
lattice is sliced along some half plane (the Volterra cut), and then deformed so that the lattice 
^ remains almost perfect away from the edge of that half plane (the dislocation line) [T5l [36] . 
Since dislocations are the principal carriers of plastic flow in crystals [23| ETJ [3l], they are 
n among the most widely studied objects of materials science. They have been investigated 

<^ analytically as points in an elastic continuum [T5| [TJ [301 [131 [31 |6], as crystal defects using 
molecular simulation techniques [HI [321 ISI], or through a variety of intermediate models such 
as phase field descriptions [T^ [211 [H]. We refer to [TTl [T^ ITE\ [B] for introductions to these 
various models. 

In the present work, we focus on the atomistic structure of dislocations. Precisely, we shall 
demonstrate that they can be understood as global minima in a variational problem. To the 
j> best of our knowledge, our results are the first that establish the existence of dislocations as 
O stable equilibrium configurations of an atomistic energy. 

^^ Our work is motivated by ongoing efforts to develop multi-scale models of dislocations such 

CN| as dislocation dynamics |6], and far-field coarse-graining techniques such as quasicontinuum 
■^ and related methods J221[3I1I33]- Our results contribute to a precise qualitative understanding 
O of the atomistic structure of dislocations, which can be used to inform the formulation and 
,__! analysis of such multi-scale schemes |^5l I2U] . 

IJ Further, our work is inspired by a recent effort to place the theory of dislocations on a rigorous 

• '-j mathematical foundation, in particular clarifying the connections between the various models 

rS mentioned above [ 211 1371 I2H1 [l2j. We outline only a small fraction of the contributions here, 

d most closely related to our own work. We believe that, to some extent, our results help to 

overcome simplifying assumptions made in many of these works. 

Possibly the most complete analysis of a static model of dislocations is provided in [T2l [TH 
[S]. This series of papers studies a continuum phase field model for dislocations in a periodic 
environment of pinning sites first described in [19j. The authors obtain a variety of scaling 
regimes, depending upon the number of the pinning sites relative to their size in terms of the 
interatomic spacing. 
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A mathematically consistent description of dislocations in an atomistic setting, using the 
language of algebraic topology is given in |1]. The concepts and language developed therein 
are used to derive models of the elastic energy of a dislocation configuration. A rigorous 
discrete-to-continuum passage within this framework is established in |28j using the language 
of F-convergence. Related works analysing dislocations and other similar defects in discrete 
systems are [U [2] . 

The works cited above concerning discrete and semidiscrete models of screw dislocations are 
primarily concerned with asymptotics of the elastic stored energy, given a number of prescribed 
dislocation cores. The creation or destruction of additional cores, for example via the intro- 
duction of dipoles, is either forbidden or explicitly tracked in the energy functional through a 
term accounting for a positive core energy. 

In contrast, our model allows for the creation and destruction of dislocation dipoles without 
any such penalty, accounting only for the stored elastic energy. More precisely, we consider a 
static atomistic model for screw dislocations in a similar vein to |28], and show that uncon- 
strained stable equilibrium configurations containing dislocations exist. Although we do not 
pursue this in the present work, we observe that the analytic properties of the equilibria we 
obtain should allow for a natural extension of the results in [28j to our unconstrained model. 

1.1. Outline. We consider anti-plane displacements of a BCC crystal, in the direction of a 
screw dislocation line (that is, parallel to the Burgers vector). For two displacements, y,y we 
consider the energy difference 

E{y;y) ■.= J2[^iDyb) - ^P{Dyb)], 

beB 

where ^B is a set of pairs of interacting (lines of) atoms, Dyi, denotes a finite difference, and 
ip is a potential describing this interaction. The potential ip is 1-periodic, where 1 is the 
atomic spacing, mimicking the behaviour of realistic pair interaction potentials. In particular, 
if y{^) — y{^) G Z for all lattice sites ^, then E{y; y) = 0. This invariance of the energy is the 
primary source of analytical difficulties. 

We call a deformation y which minimises E{y + u; y) amongst all finite energy perturbations 
u a globally stable equilibrium. 

Building on [4j and the differential displacement maps ffist employed in [35], in ^2.2 we 
present a method by which we can identify dislocation cores and assign them a corresponding 
Burgers vector. 



Our main result. Theorem 3^, states that there exists a globally stable equilibrium containing 
one geometrically necessary dislocation core. 

The proof of this result is developed throughout the remainder of the paper. In Qwe show 
that S{u) := E{y + u;y), where y is the linearised elasticity displacement field, is well-defined 



on a discrete if ^-function space describing finite energy states. In Theorem 4.5 we reformulate 



and refine Theorem 3.2 as a variational problem, by stating the existence of a global minimiser 



of S in that space, ^is then devoted to the proof of Theorem |4. 5 

Further, in ^we discuss briefly how to expoit the global stability result of Theorem 4^ to 
construct locally stable conflgurations containing flnitely many dislocation cores or conflgura- 
tions with dislocations in a domain with boundary. 

2. Geometric and Topological Preliminaries 

2.1. Anti-plane displacements of the BCC lattice. Although our analysis can be applied 
in other situations, it will be notationally convenient and physically relevant to restrict our 
attention to the body-centred cubic (BCC) lattice, which may be deflned by 

VW9 ^A/^ 
£ := BZ^ where B := s ^^ 

1/3 -1/3 1 
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and s > is a scaling factor that we leave undefined for now. 

To define anti-plane displacements, we fix the lattice vector v := [0,0, s] and define the 
projection 



which 'flattens' the Bravais lattice C onto the lattice plane with normal z/; see Figure 1(a) , It is 
straightforward to check that the set Ilu{^) is a two-dimensional triangular lattice embedded in 
M^, with lattice constant s a/8/9; we will choose s = a/9/8 so that the planar lattice constant 
is 1. We further shift the origin so that the projected set may be identified with 

A:= (^,^)^ + [ai,a2]-Z^ where ai = (l,0)^anda2= {l,^f- 

An anti-plane displacement in the direction u (or, simply, displacement) , is a map y : A — )■ 
M. The set of all displacements is denoted hj W. A displacement y gives rise to a lattice 
deformation F : £ — > M'^, 

Yiv) ■=V + y{ih "ff + n.r/)z/, r]eC. 

Let y, y be displacements and Y, Y the associated deformations. We say that y, y are equiv- 
alent if Y{C) = Y{C) (i.e., they describe the same atomistic configurations). It is easy to see 
that y, y are equivalent if and only if {y — y)(A) C Z. 

2.2. Bonds and bond lengths. Each ^ G A has six nearest neighbours, ^ + a^, i = 1, . . . ,6, 
where Oi = (1, 0)"^ and a^ = Rg~^ai where Rg denotes a rotation through angle vr/S. At a point 
^ G A, we define the set of outward-pointing nearest neighbour bonds 

7^5 := { (e, e + «.) h = 1, • • • , 6} = { (e, ^) h e A, le - r^l = 1}, 

and furthermore define the set of all bonds to be the union 

For any bond b = {C,,C, + ctj), we denote the reverse bond by —6 := (^ -|- Oj, ^). 
For b = [^,7]) & B we define the difference operator 

Dyb ■=y{r]) -y{0- 

Moreover, we set Dy := {Dyb)beB- We also note that Dy^h = —Dyb- 

With this notation, we can now define two important discrete function spaces: fixing a 
reference lattice point ^o = (0, 3 )"^, 

Wo := {v E W \ f (^0) = and snpp{Dv) is bounded}, and 

#i'2 .= 1^ g ^ I ^(^^) = and Dv G f{B)}. 

It is shown in pH Prop. 9] that W^'"^ is a Hilbert space and #0 C W^''^ is dense. 

We now introduce a crucial concept required to define the notation of dislocation. We denote 
the set of bond length 1 -forms by 

[Dy] := {a:B^ [-1/2, 1/2] | a.^ = "fe and Dyf, - ab E Z for all b E B}. (2.1) 

We note that, if Dyb ^ | + Z for all b E B, then a E [Dy] is unique, but in general there is 
ambiguity in the definition of a. This non-uniqueness is an issue which we will return to in 



^2.5 



The motivation behind this definition is that ab defines the 'shortest bond length' between 
the two lines of nuclei represented by the 2D lattice sites C,,C,', where b = {C,,^'), in that 



mm 



\y{r]) -y{v')\ = \/l + al, 
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(a) Part of £, showing the BCC unit cell, v and the (b) An illustration of the definition of a\, and its 

plane perpendicular to it. relationship to the shortest distance between atoms. 

Figure 1. Illustrations of the lattice geometry. 



see also Figure |l(b) 



The importance of the concept of bond length stems from the fact that, due to the invariance 
of the lattice under adding integer shifts to a displacement, the energy of the lattice can only 
depend on «(,, but not on Dyi, directly. 

2.3. The lattice complex. In this section, we review some terminology of discrete algebraic 
topology which is convenient for our analysis. We follow the language described in |1], where 
more details and applications to the study of dislocations can be found. 

Repeating the definitions of A, i3, we define a lattice complex as in j3l §2.3.3], with 



A:={eGM2|een.£ + (i,fr}, 

S:={(e,C)GA'|l^-Cl = l}, and 

denoting, respectively, the sets of 0-cells, 1-cells and 2-cells of the lattice complex respectively 
(see Figure [2] for an illustration). From now on, we will not explicitly use the terms p-cell, 
]?-chain and p-cochain as defined in [H §2.2], preferring instead the more evocative terminology 
'lattice points' for elements of A, 'bonds' for elements of i3, and 'cells' for elements of C We 
note the additive structure that may be defined on these objects, and write a G A to mean 
that a is an elementary p-cell contained in the p-chain A. We also frequently use the boundary 
operator 9, which maps p-chains to their boundaries, assigning orientations in the usual way. 
We then follow |H §3] in defining p-forms and integration on the lattice, writing 

eeC/ 



u 



where [/ is a p-chain, e are p-cells, and F is a p-forni (i.e. a real- valued function on p-cells). 
We note that this definition is linear in F and f/, in the sense that 

/ AF + G = A5^F(e) + A5^F(e) + 5^G(e) + 5^G(e), 

JV+V .^rr _^T/ .^TJ „^T/ 



eeu 



X F+ G + \ F+ G, 



V 
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I • Q 

■Xl-ccll ^ /*• 

/l-cell 






/ 2-cell ^ 

^ »-■■ -o 



Figure 2. An illustration of 0-, 1- and 2-cells in the triangular lattice. The 
arrows show the boundary of the positively-oriented 2-cell. Note that orientation 
only makes sense for 1- and 2-cells. 



for any A G M, p-chains U and V, and p-forms F and G. We remark here that 'bond length 



1-forms' a as defined in ^2.2 are true 1-forms in the sense defined in jH §3.1]. 



Finally, we define paths in the natural way as 1-chains 

L 

r:=$^(a,a+i), (2.2) 

fc=i 

where {C,k, C,k+i) G B for each k, and we denote the length of a path F by |F| := L. 

2.4. Measures of Lattice Distance. Since we will make use of more than simply the algebraic 
structure that a lattice complex entails, we will occasionally abuse the notation given above by 
identifying bonds and cells with their closed convex hulls; that is, we write 

X & b = {^,0 to mean x G conv{^, (}, and 

X G C = (^, C, v) to mean x G conv{^, (, r]}, 

where conv(i7) denotes the closed convex hull of a set i7 C M^; it will be clear from the context 
whether we are referring to spatial points or subchains. Since we frequently refer to them, we 
define x'" to be the barycentre of a cell C, and Cq the cell for which x'~^° = 0. 

Using this form of the notation, we define the distance from each kind of cell to the origin as 



d, 


= ie, 


db 


= inf X , 

xeb 


dc 


= inf x\, 
xec ' 



which corresponds to the usual notion of distance between sets in Euclidean space. 

The second notion of distance we will use is the graph theoretic notion. Since A can be 
identified with a planar graph with edges 6 G -B, we can further identify cells with nodes in 
the dual graph, and bonds as edges in this graph (see fU', §4.6]). This allows us to define the 
hopping distance, hop2(C, C) as the length of the shortest path in the dual graph between the 
cells C, C G C, as in [9, §1.3]. We note that since the dual graph is connected, this distance is 
always finite, and we have the following 'triangle inequality' for any dual lattice points A, B 
and C: 

hop2(A,C) <hop2(A5)+hop2(S,C). (2.3) 
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Figure 3. Two examples of bond length l-forms corresponding to the same 
deformation. The numbers are the value of the 1-form on the relevant bond, 
and arrows indicate the bond direction in which it is positive. Note the number 
and positions of the dislocation cores present change, but the sum of the Burgers 
vectors does not. 



2.5. The Burgers vector. We now define the notion of Burgers vector we use in our model, 
which is a fundamental geometric concept describing the nature of a dislocation |15| . 



We call a path T = X]fc=i(Cfc,CA:-i-i) as defined in (2.2) a Burgers loop (or, simply, loop) if 
^1 = C,L+i (or equivalently diT = 0; this implies F = ^2^ for some sum of cells A, since the 
lattice complex is perfect j5 §2.2.1, Axiom (A4)]). If F is a loop, y G W, and a G [Dy] an 
associated bondlength 1-form, then 



f a = J2Dyb + J2^ab- Dyt) = + N E Z. 



feer ber 

We call the integer A^ the Burgers vector of the bond length 1-form a around the loop F. 

Definition 1 (Dislocation Core). A dislocation core of a bond length 1-form a is a positively 
oriented 2-cell C such that Jq^jOl ^ 0. 

We will refer to cores as being 'contained in' a. For future reference, we remark that 

\dC\ _ 3 
~ 2' 



a 



dc 



< 



(2.4) 



that is, the Burgers vector around a single 2-cell can only be —1,0 or 1, and hence we define 
the sets of dislocation cores 



C+H := I 



C eC 



C positively oriented, / a 

dC 



C-[a\ -.= 10 eC 
C^[a] ■=C^[a]UC-[a]. 



C positively oriented, / a 

dC 



+1 
-1 



Remark 2.1. It is interesting to note that if a, a' G [Dy], then they need not have the same 
number of cores; see Figure [3] for an illustration of this fact. 

The only point at which this ambiguity is an issue is if a has C, C G C^ [a] which are 
adjacent. In that case, it may be checked that the h G dC such that — 6 G dC must have 
ah G { — 1/2,0, 1/2}. In the case where a^ = ±1/2, redefining at = =Fl/2 removes these cores, 
and a remains a bond length 1-form in [Dy], so we will always assume that minimising sequences 
have Ob = for any bond b shared by 2 adjacent cores. D 
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The net Burgers vector is obtained by summing the signs of the cores, or equivalently, by com- 
puting the Burgers vector on a sufficiently large loop enclosing all cores. Since a G [Dy] is not 
necessarily unique for a given y, we ensure that such a concept can be defined unambiguously. 

For our purposes it will be enough to consider displacements with some prescribed far-field 
behaviour. 

Proposition 2.2. Let y E W and a G [Dy] such that a^ — )■ as (i^ — )■ oo. Then, for any 
a' G [Dy], C^[<y'] is finite and 

E / " = E /«=/«' (2-5) 

where T is any loop that encloses all cores in a. 
Proof. If ttf, — > as (i;, — i- oo then 

/ a — )■ as dc —> oo. 

JdC 

Since J^^ a G Z it follows that J^^ a = for dc sufficiently large, and hence the number of 
dislocation cores present in a is finite. 

Moreover, since a^ G (—1/2,1/2) for di, sufficiently large, it follows that Of, = a'j, for all 
a' G [Dy] and d^ sufficiently large. In particular, C^[a;'] is also finite. 



To prove (2.5), let F be a loop that encloses all the cores in a for which F = dA. Then 



a = y^ / a = y^ I a 

^ C&A "^^<^ C£C± ^^ 



Taking F such that Dyi, = ab = a^ for all 6 G F we obtain the ffist identity in (2.5) as well. D 
We can now formally define the net Burgers vector. 

Definition 2 (Net Burgers Vector). Let y E W such that ab ^ as dh ^ oo for some 
a G [Dy] . Then we define the net Burgers vector of y to be 



B[y]:= E / 



a, 



cec± fal""^^ 



for an arbitrary a G [Dy] . 

The quantity B[y] can be experimentally observed from outside the system, by determining 
the strain at 'infinity'. For example, if i?[y] = 1, then this tells the observer that there must 
be at least one dislocation in the system, but nothing about the total number. 



3. Main Result 
In this section we present the main result of the paper with accompanying assumptions. 

3.1. Energy difference functional. Before we can state the main result, we introduce an- 
other key concept that we employ in its formulation and proof: the energy difference functional. 
We assume that lattice sites (corresponding to lines of atoms in the BCC crystal) interact via 
a nearest-neighbour pair potential ip G C{M.) fl C'*(]R \{Z + 1/2)), which satisfies the following 
properties: 
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(ipl) ip is 1-periodic; 

{■ip2) ip and ip{^ + ■) are even; 

(■j/^S) iplr) = if and only if r G Z; 

(V^4) ^"(0) = /i > 0. 

(V'5) ^{x) > ltp"{0)x'^ for all x G [-^, ^]. 

Remark 3.1. Assumptions {ip\) — (V'4) are very general, and are natural in the physical 
context: {ipl) and {ip2) encode lattice symmetries, while (t/^S) and {ip-i) state that the system 
has a stable crystalline ground state. 

The only "technical" assumption is (V^5). The reason for this assumption will become apparent 
in §5.1[ where we use it to establish an a priori bound on the number of dislocation dipoles in 
finite energy configurations. We believe that (^/'S) can be replaced with weaker variants, but 
cannot be removed altogether. 

We remark that the requirement that ip G C'^(]R \ (Z + 1/2)) can be relaxed further by 
modifying the proofs we give below, but since this adds little at the expense of readability, we 
omit such arguments here. The prototypical example of a potential satisfying {ipl) — ('?/'5) is 
^(r) = ^'""(r) := |dist(r, Zf. U 

For two displacements y,y E W we define the energy difference functional, formally for the 
moment, as 

E{y-y):=Y,[i^{Dy,)-ij{Dy,)\. 

beB 

For example, iiy — y G Wq, then E{y; y) is clearly well-defined since the sum is effectively finite. 
For arbitrary displacements y, y, E{y]y) need not be well-defined. However, we will show in 



§4.1| that E can, under certain conditions, be extended by continuity to relative displacements 
y-y e W^'^. 

Using the terminology of energy differences, we can define what we mean by a stable equi- 
librium displacement. Intuitively, the definition entails that finite energy perturbations cannot 
lower the energy. 

Definition 3 (Stable Equilibrium). A displacement y E W is a locally stable equilibrium 
if there exists e > such that E{y + u;y) > for all u eWq with ||-Dm||2 < £■ 

A displacement y eW is a globally stable equilibrium if E{y + u;y) > for all u G #o- 

3.2. Statement of the main result. Recalling Definitions [2] and |3] the existence of a screw 
dislocation can be formulated as follows. 

Theorem 3.2 (Existence of a geometrically necessary dislocation). There exists a 
globally stable equilibrium displacement y eW with net Burgers vector B[y] = 1. 

The notion of global stability described in Definition [3] is equivalent to the statement that 
a displacement is stable if any finite energy perturbation increases the energy of the system. 
Describing a dislocation configuration as the minimiser of an energy difference functional gives 
us access to the Direct Method of the Calculus of Variations. 

We refer to this result as the existence of a 'geometrically necessary' dislocation since we do 
not prescribe the absolute number of dislocation cores, but only the net Burgers vector. 



Outline of the proof of Theorem 3.2 



We define a reference configuration y(^) = ^ arctan(|^) (the continuum linear elasticity 
solution for a dislocation), with the aim to minimise the energy difference functional 
£{u) := E{y + u; y) over a suitable class of functions u. 

We show that this functional, initially defined over >^, can be continuously extended 
to a functional over W^'"^. 
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(2) In order to use the Direct Method to estabhsh the existence of a minimiser to £, the 
crucial step is to obtain a global lower bound on the energy. This is the main step in 
the proof, and requires careful geometric estimates based on the number of dislocation 
cores and the distance between them. We shall prove that £{u) > \\fi\\\ — 1, where (3 



can be thought of as belonging to [Du] (however, see (5.1) for the precise definition). 

(3) This lower bound guarantees in particular that the number of dislocation cores is 
bounded along a minimising sequence as well as weak compactness of a minimising 
sequence m". 

(4) The final step is to ensure that limw" has non-zero net-Burgers vector. This need not 
be the case since weak convergence of Du^ allows for energy to be translated to infinity. 
In our present context it is possible, by introducing a dislocation dipole, to effectively 
translate the geometrically necessary core to infinity, and thus obtain a limiting dis- 
placement with zero net Burgers vector. We shift the minimising sequence and employ 
a concentration compactness argument to prevent this. 

D 



3.3. Locally stable equilibria. Theorem |3.2| establishes the existence of a configuration y = 
y + u, which is a globally stable equilibrium configuration for a single screw dislocation in an 
infinite lattice. From this starting point, it is possible to construct more general locally stable 
equilibrium configurations. The idea is (1) to superimpose copies of y and define 



m--=J2'^yi^ -'"''')■■ 



where Cj G C are cores in z and Sj G {±1} the Burgers' vectors of these cores; (2) to show that 
z is an approximate equilibrium when the cores Cj are sufficiently far from one another; and 
(3) to apply the inverse function theorem to establish the existence of an equilibrium z close 
to z. 

Here, we only state two results that we obtain by this strategy, but refer for their proofs 
to jTTj, where we present them in a more general context. 



In step (3) of the strategy outlined above we require a discrete ellipticity condition (3.1), 
which can be established rigorously, for example, for a piecewise quadratic potential. 

Lemma 3.3 (Discrete ellipticity). Let ^(r) := ipunir) := |dist(r, Z)^ and let y = y + u, 

u G W^'"^ , be a locally stable equilibrium configuration. Then Dyh G M \ (^ + Z) for all b E B, 
and hence 



J2^"{Dyb)Dvl>X\\Dv\\l Vt; G #''^ (3.1) 



b&B 



Proof. Suppose there exists a bond b = (^, C) such that Dyt E ^ + Zi. Let Zt{ri) := y + tSr^^^, 
then a direct calculation shows that E{zt;y) < for some sufficiently small t (either positive 
or negative). D 



The two results we state in the following admit general ip, but require (3.1 ) as an assumption: 



(STAB): There exists a locally stable equilibrium y = y + u, u E W^'"^, satisfying the 



ellipticity condition (3.1 ). Moreover, let A be a finite union of cells such that C^ [a] C A, 
for any a E [Dy]. 

Our first local stability result states that any configuration of dislocations is stable pro- 
vided that the cores are sufficiently separated. In particular, it shows that there exist stable 
configurations with arbitrary net Burgers. 

Corollary 3.4 (Finitely many cores). Suppose that (STAB) holds. 
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Let Cj & C,j = 1, . . . J, J ^ N, be a finite collection of cells, and let Sj G {±1}- There is a 
minimal separation distance Lq > such that, if miiii^j \x^^ — x'"'\ > Lq, then there exists a 
locally stable configuration z & W such that, for any a G [Dz], 

J 
C^[a] C [j {x^^ + A) and 

i=i 



d{x^:i+A) 



a = Sj, forj = l,...,J. 



In particular, B[z] = Ylj=i^ 



]■ 



Our second local stability result states that dislocations are stable provided they are suffi- 
ciently distant from any domain boundary. To state this result, let fi := {.^ G A | ^2 < 0} be a 
discrete half space and let B^ := {h = {^, Q E B \ S,-, C, (z Vt} be the corresponding set of bonds. 

Corollary 3.5 (Domain with boundary). Suppose that (STAB) holds. 

Let C E C such that L := — (x'")2 > 0. If L is sufficiently large, then there exists a locally 
stable half-space configuration z : Q -^ M. containing a dislocation. That is, for any a G [Dz], 

(1) C^[a] Cx^ + A, 

(2) /a(.c+A)« = l. ^nd 

(3) there exists e > such that 

Y, [^Dzb + Dvb) - ^{DzS) > Vt; G #^o, ll^^lh < e. 

3.4. Regularity. The globally stable equilibrium configuration y, whose existence we estab- 
lished in Theorem 3.2 is of the form y = y + u, where y(^) = ^arctan(|^) is the continuum 
linearised elasticity solution for a screw dislocation. We refer to Qand in particular to Theorem 
14.51 for further details. 

This fact implies that only a finite amount of energy is stored in the dislocation core, and 
that, up to some fixed prescribed error tolerance, the linearised elasticity displacement field is 
accurate outside of some fixed radius. These observations give rise to new points of view on the 
concepts of dislocation core energy and core radius, which we explore in [I6j . In particular the 
core radius is an interesting concept related to the decay of the "corrector" u to the configuration 
y predicted by linear elasticity. Here, we state a regularity result proven in more general form 
in [IQJ, which precisely quantifies the rate of decay of Du. In effect, the results states that the 
decay of Du is the same as predicted by linearised elasticity. 

Proposition 3.6. Let y = y + u,u E W^'^ be a locally stable equilibrium, then there exists C 
such that 

\Dub\ < Cd^^ V6 G B. (3.2) 



Remark 3.7. One may expect, and numerical simulations confirm this, that the corrector 
u satisfies the three-fold symmetry of the lattice A with respect to its origin (recall that the 
origin lies in the barycentre of a cell). Exploiting this symmetry, one can observe that the 
decay rate is in fact \Dub\ < Cd^^. However, as soon as the symmetry is broken, for example 
by applying a small shear displacement at infinity, or by moving the core off the centre of the 



cell, the generic rate (3.2) is observed also numerically. D 
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4. Analysis of the energy difference functional 

4.1. Extension of the energy difference functional. We fix a displacement y and define the 
functional £{u) := E{y + u;y). For u & Wq this is always well-defined. If Dyt G M \ (Z + B{e)) 
for all b & B, where e > 0, then the first and second variations (in the sense of directional 
derivatives) are also well-defined, and given by 



{5£{0),v) = J2'P'iDyb)-Dvb, Ioxv^Wq, and (4.1) 

haB 

(fE{^)v, w) = ^ fiDyb) ■ DvbDwb for v,w e Wq. (4.2) 

beB 

S'^S{0) can clearly be extended by continuity to v,w ^ W^'"^, but this is less obvious for 6£{0) 
or for S itself. We first state a general result. 

Lemma 4.1. Let y & W satisfy Dyh EM.\{'L + -B(e)) for some e > and suppose that 6£{0) 
is a bounded linear functional ({6£{0),v) < C\\Dv\\2 for all v G #0^- Then, £ : Wq ^ M. is 
continuous with respect to the norm \\D ■ II2; hence, there exists a unique continuous extension 
of£ toW^'\ 

Proof. The proof of this result is analogous to the proof of Theorem 2.8 (ii) in [26]. For 
convenience we give a brief outline. 

For M G y^ it is easy to see that 

£{u) = Y, bf^iDyb + Dub) - ^{Dyb) - ^'{Dyb)Dub\ +Y,^'iDyb)Dub. 
beB beB 

Since we assume that 6£{0) is a bounded functional, the second term on the right-hand side is 
continuous. Using the fact that ||Dw||oo < ||-Dw||2, the smoothness oiip, and the fact that each 
summand in the first group is effectively quadratic in Dub, it is easy to show that the second 
term on the right-hand side is continous as well. D 

For future reference, we now derive a simple condition on y under which 6£{0) is a bounded 



functional. Applying summation by parts to (4.1) we obtain 



{6£{0),v) = J2fiO-m, where /(O := $^ ^'(/^yb(0) (4-3) 

is the force acting on atom ^ under the displacement y. The following result states that, if y is 
sufficiently close to equilibrium in the far-field, then 6£{0) is a bounded linear functional. 

Lemma 4.2. Suppose that a displacement y has associated forces /(^) satisfying the bound 
1/(01 < Ci{l + 1^1)-* for some t > 2, then {6£{0),v) < C2\\Dv\\2 for alive Wq. 



Proof. Proposition 12 in ^24j immediately implies that 

IIp^^-^II <C\\Dv\\2, (4.4) 

II log(l?l+2) lloo — II IM' V / 

for some constant C > 0. (This inequality is essentially a consequence of the embedding 

IIi'IIbmo < C\\Vv\\l2 for v G C\R^).) 
We can therefore estimate 



|(5£(0),t;)|<5^|/(Olb(OI<||log(|ei + 2)/||, 



log(|5|+2)lloo- 
«6A 



The assumption |/(0| < C|^|-* with t > 2 implies that || log(|0 + 2)/||i is finite. D 
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Cut 



Figure 4. An illustration of the coordinate system and reference configuration chosen. 

4.2. The reference displacement. We now specify the reference displacement y used in the 
definition of S in §4.1[ It is best to think of y as prescribing a far-field boundary condition 
y{0 ~ yiO ^s 1^1 ~^ ^^- We wish to choose y in such a way that it enforces a geometrically 
necessary dislocation, and at the same time satisfies the condition of Lemma |4.2[ 



A natural choice is the dislocation displacement field from linear elasticity theory. Since it is 
instructive (though not essential to our proofs) we give a brief motivation of this construction. 
In the far-field, we expect that continuum linearized elasticity theory is a good approximation 
to the atomistic equilibrium condition 6S{0) = 0. This can be formalized by first deriving the 
Cauchy-Born approximation and then linearising it. Due to the hexagonal symmetry of A one 
finds that the linearised continuum approximation is simply Laplace's equation, Ay{x) = 0. 

Hence, following Section 3-2 in Hirth & Lothe [15], we define y as follows: 

y{x) := ^ arg(x) = ^ arctan (f^) , (4.5) 

where we identify x G M^ with the point Xi + 1X2 € C, and the branch cut is taken along the 
positive ^i-axis, as shown in Figure |4j 

The gradient (away from the branch cut) is given by 



v27rr2' 27rr2. 

where r := |x|. This function can be extended to a function in C°^(]R^ \ {0}), which we take as 
the definition of Vy from now on. Moreover, we can check that indeed A^(x) = div(V^(x)) = 0, 
in the pointwise sense, for x 7^ 0. 

Let a = {ab)b&B be a bond-length l-forni associated with y; we claim this is unique, and the 
following lemma provides a convenient formula for a^ in terms of Vy. 

Lemma 4.3. Let a G [Dy] then for any bond b = (^, ^ + ctj) G B, we have 

&b= I S/y{i + tai) -Oidt. (4.7) 

Jo 

Proof. By definition, Vy is independent of the choice of branch cut. Moreover, if the branch 
cut is chosen differently, then the displacement at each site is only changed by an integer, which 
means a does not change; hence, a is also independent of the branch cut. 

Now fix 6 = {^,^ + ai) G B. Since the origin lies at the centre of a cell we can redefine y with 
a branch cut that does not intersect b. The Fundamental Theorem of Calculus gives 

/ Vy (e + tOi) ■ a, dt = y{i + a,) - y(0 = ^ ( arg(e + a.) - arg(O) , (4.8) 

Jo 

and since we have assumed that 6 is a nearest neighbour bond, it has length L The term on 

the right hand side is l/27r times the angle formed by the points ^, and ^ + Oj, which is 
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maximised by making ^ and ^ + Oj as close to the origin as possible — that is, when ^ and 
^ + Oj are on the boundary of Cq. It follows that the angle can be no larger than ^, and hence 
c^b = vi^ + CLi) "ViO ^ ["!' |]- This implies that a is unique, since Dyb ^ ±| for all b E B. D 

As an immediate corollary of Lemma [4. 3| we obtain the following bound on d;,: 

1 



\ab\ < 



27idh 



ybeB. 



(4.9) 



We conclude the analysis of y by showing that it satisfies the conditions of Lemma 4.2 



Lemma 4.4. Let y be defined by (4.5), and let /(^),^ G A, be the associated forces (see 



(4.3 )j, then a G [Dy] satisfies a^ G [—1/3, 1/3] and 

1/(01 <ier' v^ga. 



In particular y satisfies all conditions of Lemma 4-2 



Proof. Recall from (4.3) that 



benio ben(e.) 

Taylor expanding ip'f^ to third order, using the fact that ^ip'^O) = 4''"{0) = (since ip is even 
about 0), gives 



/(0= E [r{0)&b+l^^'\s,){a,) 
benio 

for some Sb G conv{0, df,}. Applying ( |4.9[ ) we obtain 

/(0= E r{0)ab + O{d^'), 
benio 



(4.10) 



We now inspect the sum on the right-hand side of (4.10) in more detail. Applying (4.7) we 
rewrite this sum as 



6 „i 

E V'"(0)d6 = E ^"(0) / ^^(^ + ^«0 ■ «' d^- 

Taylor expanding Vy(^ + taj) and using the fact that |V'*y(a;)| < \x\~^, we obtain 

/ Vy{^ + tai)■a^dt = Vm■a^ + lV^mK(^i] + |V'y(0[a„ a^, a^] + 0(|er'). 
Summing over i = 1, . . . ,6 the first and third terms cancel since a^+s = — a^, hence we obtain 

E riO)ab = li2^^'^'y(Oa^ + 0{\^\-'). (4.11) 



fce7^(0 



i=l 



We now observe that 



lY.^IV'y{Oa^ = -|Am = 0. 



i=l 



Inserting the last identity into (4.11 ) and combining the resulting estimate with (4.10 ) we obtain 
the stated estimate on |/(^)|. D 
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4.3. The variational problem in W^'"^. Combining Lemma 4.4 with Lemma 4.1 and Lemma 



4.2 



we deduce that £{u) := E{y + u\ y) is a well-defined and continuous functional on ^^'^, 



where y is the reference configuration defined in (4.5). It will later be convenient to recall from 



the proof of Lemma |4.1| that the explicit definition of the extension is 

£{u) = J2 [i^ic^b + Dub) - i^i&h) - i^'iab)Dub] + {5£{0),u), (4.12) 

beB 

In the next section, § |5} we will prove the following result: 
Theorem 4.5. There exists u E W^''^ such that S{u) < £{v) for all v E W^''^. 



As an immediate corollary we can now prove Theorem |3.2 



Proof of Theorem \3.^ Let y := y + u, where u is a minimizer of S in #'^'^. Since Du E (?{B) it 
follows that \Dui,\ — )■ uniformly as db — ?■ oo. Using also the fact that |db| — )■ uniformly (cf. 



(4.9)), we conclude that a^ = CH) + Dub + Zb, where Zb is a compactly supported, integer-valued 



1-form. From the definition of the net Burgers vector and from (2.5 ), it now follows immediately 
that B[y] = B[y] = 1. Moreover, minimality of u implies that y is a globally stable equilibrium 
in the sense of Definition |3l D 



Theorem 4.5 is interesting in its own right: it shows that atomistic configurations containing 



dislocations can be obtained as global minimizers of a variational problem formulated over A. 



This is particularly useful for further study (e.g., of regularity; cf. § 3.4) of dislocations in this 
model. 

We also remark that any local minimizer m of £^ in W^'"^ would give rise to a locally stable 
equilibrium with net Burgers vector B{y + u) = 1. The advantage of local minimisers is that 
they can be computed numerically. 

5. Proof of Theorem 14.51 

As currently formulated, it is not obvious that the energy S is bounded below, and it is even 
less clear whether S is coercive in a sense which would allow us to invoke the Direct Method. 
This is due in large part to the fact that the reference configuration is nonlinear and ijj is 
periodic, so the integrand has infinitely many energy wells. 

The periodicity of ip allows the creation of dislocation dipoles 'cheaply'. If dipoles are well- 
separated, then each dipole gives a positive contribution to the energy which is proportional to 
the logarithm of the dipole length (the separation distance between the two cores of the dipole). 
However, for generic configurations of dipoles the sign of the energy contribution is difficult to 
determine, since it depends strongly upon the relative orientations of the dipoles. In essence, 
this is a geometric nonlinearity of the system, and most of the effort expended in what follows 
will be to control the number of dipoles that can form. 

From a technical point of view the issue arises as follows: in §2.2 we decomposed Dy = a + w, 



a E [Dy] , since the energy of the displacement y only depends on a due to the periodicity of 
the potential tp. Consequently, if we have a sequence u"^ with S{u^) uniformly bounded, then 
this will bound only ||/3'^||2 for P"' £ [-Du"], and not ||Du"||2. In particular, generic minimising 
sequences cannot be weakly compact. 



By exploiting the vertical shift invariance (^5.2) and the horizontal translation invariance 



5.3| of the energy £, we will construct a weakly compact minimising sequence. Having made 



this special choice of minimising sequence, we use a profile decomposition in ^5^ We show that 
each profile obtained in this way has net Burgers vector zero, leading to the conclusion that 
the net Burgers vector of the limit remains 1, and proving existence of a minimiser with the 
properties required. 
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5.1. An elementary lower bound. Our eventual goal is to establish a coercivity result for 
£. We begin with an elementary lower bound that will motivate subsequent constructions. 

Let y = y + u,u & W^'"^, be a trial displacement, a G [Dy], and recall that a = [Dy] is 
unique. Since u G W^'"^, and hence Dub — > as (i;, — )■ oo, it follows that y has a well-defined 
net Burgers vector in the sense of Definition [21 and B[y] = B[y] = 1 . 

Let 

(3 := a — a; (5.1) 

this 1-form satisfies the property that 



/3 7^ if and only if a j^ a, 

dC JdC JdC 

that is, dislocation cores present in /3 are those that are introduced by the addition of u to y. 

Remark 5.1. We note that /3ft does not necessarily belong to [—1/2, 1/2], and hence is not 
a bond length 1-form, so the definitions of §2.5 do not strictly apply; however, it remains a 
1-form in the sense of [H §3.1]. As L^^/? G {0, 1, —1} for all C G C, we shall therefore slightly 
abuse our notation and refer to dislocation cores in (i as the cells C ^ C for which 



/3 = ±1. 

'dC 

We also define C+[/3], C~[f3] and C^[/3] in the obvious way. D 

Next, we define z : B ^ Z via 

Du = /3 + z, (5.2) 

which is compactly supported since P,Du G £'^{B). We shall see in ^5.2 that the support of z 
can be thought of as a union of branch cuts connecting dislocation dipoles. 
With this notation, we obtain the following result. 



Lemma 5.2. For any u G W ' with Du = (3 + z as in (5.2) and for any e > 0, we have 

£{u) > (|^"(0) - e) m\% - J2 fi^b)zt + {5£i0), u) - C,, (5.3) 

fees 
where C^ > is a constant that is independent of u. 



Proof. This estimate arises from the expression (4.12); using the periodicity of the potential ijj., 
we can write 

^(^) = I] (^(«fe + Pb) - ^("'') - ^'(afe)A) - I] ^\(yh)Zb + (5f (0), m). 

Define the function 

^{t^s)-^{t)-i,\t)s 
9{s.t):={ J. ^^0' 



2 



i^"(t) s = 0. 



By assumption ('?/;5) in ^3.1, (7(5,0) > -^i\)"{^) for any |s| < 1/2. 

Since g is uniformly continuous on [—1/2, 1/2] x [— r, r] for some r > 0, it follows that for 
each e > there exists (5(e) > such that 

c/(s, t) > ^V^"(0) - e for |s| < ^ + S{e) and |t| < (5(e). 

Next, we note that (4.9) implies 

1 1,^,1 ^ , 1 1 

\ab\<^;—r and \f3b\ = \ab - ab\ < - + -—r- 
zirdb 2 Indb 

Hence there exists Rq > such that, for db > Rq, 

9{/3b,ab)>W'{0)-e, 
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which can equivalently be stated as 

^Pi&b + 13,) - ^{ak) - ^\ab)Pb > (|^"(0) - e) \Pb\^ for 4 > Ro. (5.4) 

It may be checked that 

i^{b\d,<Ro}<Rl 
and since ip, ip' and (5 are uniformly bounded, it therefore follows that 

beB beB 

We can think of ||/3||2 as estimating elastic stored energy. In the following sections we will 
establish several results on z = Du — /3, which will eventually allow us to bound the remaining 



terms {S£{0),u) and X^b'^'l^b)-^^ i^ (5.3). 



5.2. Dipoles &; Branchcuts. Let y = y + u, u G W ' , be a trial displacement, a G [Dy], and 



let /3, z be defined by (5.2). While a and hence j3 are uniquely defined (except in borderline 
cases when a, € {±1/2}), one can exploit the vertical shift invariance of the lattice (encoded 
in assumption (ipl), periodicity of ip) to construct equivalent displacements u G W^''^, 

u:=u + U (5.5) 

where U : K ^"L and t/ G #o; and hence modify the z component. 

If we let ^ := y+u, then clearly, a G [Dy\ and this leads to the same definition of /3. Crucially, 
though, Du — /3 7^ Du — /3. We can therefore ask how to choose U in an "optimal" way. It 
turns out that minimizing the total length of the branch cuts is a useful choice, which amounts 
to minimizing \\Du + DU — /3||i = \\Z + D[/||i. Since z has compact support, a minimizer 
clearly exists, but it need not be unique; see Figure [5j We may therefore assume, without loss 
of generality, that u satisfies the discrete minimal connection property (DMCP) 

||L>n-/3||i = ll^lli = min \\Du + DZ - (3\\i. (5.6) 

This minimality condition is similar to the idea of minimal connections, introduced in [5]. 



We will now establish various properties of the structure of z defined in (5.2). In particular, 
we will show that z can be decomposed into a sum ^ z"^ and that the support of each z™" is 
analogous to a branch cut for a dipole. 



Lemma 5.3. Let u G W^'"^ satisfy the DMCP (5.6) and suppose Du = P + z as in (5.2). 
Then we can write 

M 

z = J2'"'^ (5.7) 

m=l 

where M = #C"'"[/3] is the number of dipoles contained in (5 and z™' : B ^ { — 1,0,+!}, 
m = 1, . . . , M , satisfy the following properties: 

(1) zjp^ = 1 on a sequence of bonds (6j)"=o such that 

(a) dbi and dbi+i share a common 0-cell for each i = 0, . . . ,n — 1, 

(b) bo G SC^ and -bn G dC+ where C+ G C+[/3] and C^ G C-[f3]. 

(2) 2;™ = for bonds outside the set {±6j \i = 0, . . .n}. 

Proof. The result is geometrically intuitive; see Figure |5j We therefore postpone a complete 
proof to Appendix |A.l D 
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Figure 5 . A typical example of the support of a mininial z and its decomposition 
into z"^ for a given distribution of dipoles. An arrow means that z^ = 1 on the 
bond pointing in the direction of the arrow, and a double arrow means Zi, = 2. 
The grey dashed bonds on the left of the diagram show an alternative definition 
of z with the same minimal norm, and the black dashed lines show two lines 
which show that z"^ can be decomposed into 2 straight cuts. 



We will say that z"^ connects the dislocation cores C^ and C„. Moreover, we obtain the 
following corollary. 



Corollary 5.4. Each 1-form z^ in the decomposition (5.7) can he identified with a shortest 



path in the dual lattice between C^^ G C^ [a] and C^ E C [a], and further 

||^'"||l=hop2(C+,C-). 



Proof. First, we observe that, due to the DMCP (5.6) and the decomposition proven in 



Lemma 5.3 



P 1 



M 

E 

m=l 



(5.8) 



That is, if z^ 



b ; ^b 



7^ 0, then z^ 



b ; ^b 



have the same sign. 



The construction employed in the proof of Lemma 5^ identifies a sequence of bonds bi and 
cells Cm,i such that bi, — fej+i G dCm,i, ^fe" = 1- Using the natural identification of cells with 
points in the dual lattice, this implies that Cm,i are adjacent in the dual lattice, and furthermore 
that bi can be identified with edges connecting these cells; this leads to the fact that 

||^™||l>hop2(C+,C-). 

To prove the converse, we take the path in the dual lattice corresponding to z'"^, adjoin a 
shortest path between C^ and C~ in the dual lattice, and thus obtain a closed dual lattice 
path. We construct a polygonal closed path in M^ by connecting the barycentres of the cells 
along the path and define U to be the characteristic function of the bounded interior of this 
loop in M^. 

It is now straightforward to check that by defining z"^ := z"^ + DU and z := z — z"^ + z"^ 



M 
m'=l 



\z 111 + If 
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Since f/ is a compactly supported integer shift as in (5.5), the discrete minimal connection 



property (5.6) implies that ll^z^lli < ||-S™||i = hoP2(C*m7 C*m)) completing the proof. D 

Later on it will be convenient to assume that each cut z"^ is made up of at most two straight 
cuts: straight cuts are defined to be 1-forms z : B ^ { — 1,0,1} for which there exists a line 
L := {x'" + taj 1 1 G M}, where x^ is the barycentre of some C & C and Oj a nearest neighbour 
direction, such that Zf, ^ ii and only if the bond satisfies 6 fl L 7^ 0, 

clos{x eR'^\xebeB,Zb7^o} 

is a connected set, and ,2^ > either exclusively on bonds in the directions Oj+i and 0^+2 or in 
the directions ai_i and ai-2- We will say that a straight cut 'lies in the direction a/ whenever 
ttj = ±aj in the definition of the corresponding L. See the cut depicted in the centre of Figure |5] 
for a visualisation of the definition. 



In the next lemma, we show that we can always choose the decomposition (5.7) such that 



each z"^ is composed of at most 2 straight cuts. We will refer to any u as in the conclusion of 



Lemma 5.5 as satisfying the straight cuts property. 



Lemma 5.5. Let u E W^'"^, and Du = (3 + z as in (5.2). Then there exists u E W^''^ 

satisfying the DMCP (5.6) as well as Du = (3 + X]m=i ^™ where each z^ is the sum of at 
most 2 straight cuts. 

Proof. The idea is to show that we may always find a shortest path in the dual lattice between 
any pair of cells which is made up of 2 straight segments. It is intuitively clear from Figure |5] 
that this can always be done. 

A complete proof is postponed until Appendix |A. 2 D 



5.3. Shifting the Origin. Suppose that y = y + u, where u E W^'"^ satisfies the DMCP (5.6). 
For any C G C, we define the affine transformation 

c* _ / ^ + x'" if the triangle C points upwards, , . 

\ ^6(^ + 3;'") if the triangle C points downwards, ^ ' ^ 

where Re denotes the rotation through angle vr/S. Thus, F*^ maps the lattice onto the lattice, 
and Co onto the cell C. 
Next, let 

u^ ■=uoF^ + yoF^ - y. (5.10) 

It follows that y(^) + u'"{C,) = y{F^^) + u{F^^) for all lattice points ^ G A, so that there are 
corresponding bond length 1-forms a'" E [Dy + Du'"] satisfying 

a^ = a o F^. 

As before, define /3*" := a'" — a. According to these definitions, 

yiF^'O - m) ■■=l{ arg(e + x^) - arg(O) ; 

if we make this function single-valued by introducing a compact polygonal branch cut pass- 
ing through the barycentres of a shortest dual lattice path between Co and C, then it is a 
straightforward exercise to show that 

\D{yoF%-Dy,\<d^'; 
therefore u^ E W^'"^, and 

£{u'') = E{y + u'';y) 

= E{yoF^ + uoF'^;y) 

= E{y oF^ + uoF^;yoF^) +E{yo F^; y) 

= S(u), 
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noting that the the first term on the third hne is simply a resummation of S{u), and the second 
term vanishes. 

For each C G (^"""[a], we can replace u'" with u^' = u'" + U for some f/ : A — )■ Z, such that 
\\Du'^ + DU - /3'^||i is minimal, i.e. u^ satisfies the DMCP ( |5^ . We obtain that 



To summarize, we have constructed a corrector displacement u'^' 



'1,2 



with the same energy 



as M, but for which C G C"'"[a] has been shifted to the origin. Upon minimising \\Du — (3 
amongst all choices C G C+[a], we obtain the following result. 



ic\ 



Lemma 5.6. Let v G W^'"^ , then there exists u G W^'^ such that £{u 
the discrete optimal connection property holds: 
(DOCP) There exists a G [D{y + u)] such that 

\\Du-[3\\^= min min \\Du^ + DU - P^h, 



£{v) and such that 



where (3 



a 



a, u is defined by (5.10) and (3' 



mm 

C£C+[a]U:A^Z 



a o 



(5.11) 
F'^ — a, where F'^ is defined in (|5.9). 



The crucial property that we obtain from the DOCP (5.11) is a bound on the distance 
between the necessary core at Cq and all negative cores. 



Lemma 5.7. Suppose u G 



'1,2 



satisfies the DOCP (5.11) and let z = ^ 






to (5.7). Then, 



according 
(5.12) 



hop2(Co, C^) > hop2(C+, C^), /or m = 1, . . . , M, 
where we recall that z^ connects the cores C^ G C"*" [a] and C^ G C~ [a] . 

Proof. Suppose the converse for contradiction. Then there exists m and a dual lattice path 
connecting Cq to C^ which is strictly shorter than hop2(C;J^,C~). Letting F := F*^™ and 



V :- 



u 



^77 



Dv 



^oF = zoF + S = ^z"'oF + z, 



where z is the contribution coming from the branch cut in y o F — y. Consider the closed 
curve passing from FC^ = Cq to FCq along the branch cut, then along a shortest lattice path 
between FCq and FC^, and then back to FC^ along the support of z"^ o F. This is a closed 



curve, and by a similar argument to that in Corollary 5.4, we can define w & Wq as w{^) = 1 for 
^ G A inside the curve, and outside. It can then be checked that v + w has a corresponding 
z which satisfies 

||^||i = ||2;oF|| - ||2;™oF||+hop2(Co,C„) = ^ H^^li + hop2(Co,C^) < \\z\\i, 

the required contradiction. D 

As a corollary we obtain the following stronger property. 



Corollary 5.8. Suppose u G W ' satisfies the DOCP (5.11 ) and let z = X]m=i ^™ according 



to (5.7). Then, for any m G {1, . . . , M} and for any cell C G C such that zj^ ^ for some 
bedC, 

hop2(Co,C)>hop2(C+,C). (5.13) 



Proof. Lemma 5/7 states that, if C G C and C^ G C^ are connected by z™, then 

hop2(Co,C-)>hop2(C+,C-). 

It is clear that any subpath of a shortest path in a graph is also a shortest path. By the 
construction of z™, C lies on a shortest path between C^ and C~ , and therefore 

hop2(C+, C-) = hop2(C+, C) + hop2(C, C-). 
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The triangle inequality for paths (2.3) now directly implies 



hop2(Co,C)>hop2(C+,C). 



D 



5.4. Estimating {S£{0),u). In ^5.2 and g5.3| we showed that for any u G W^''^, we can find 



u G W^'^ such that S{u) = S{u), and for which the corresponding branch cuts z satisfy 
the DOCP(5.11). We are now in a position to exploit the chosen structure of z to derive 



compactness for minimising sequences. 

Our first step is to provide a stronger bound on S£{0). We have already shown that 
\{6S{0),u)\ < 11/^1x11^2, but this will not be sufficient since our estimates so far only provide 
a bound on /3, and not on Du itself. Therefore we need to estimate [(^^^(O),^)! only in terms 

of |l/5||2. 



Lemma 5.9. For each u E W^''^, let f3u '■= (3 be defined through (5.1) 
There exists a constant C > such that 



[6S{0), u) < CWPuh Vm G ^^'=^ satisfying the DOCP ( |5lT| 



We provide the proof of this fundamental estimate throughout the remainder of this section. 



Recall the definition of ^o from §2.2[ For any ^ G A, it is always possible to express the 
difference ^ — ^q as 

i-io = nai + mai+i, 

for some nearest neighbour lattice direction at and some n,m G N U {0} with n 7^ unless 
^ = ^0- We then define the path F^ to be 

n—l m— 1 

F^ := J^ (^0 + jai, ^0 + U + l)ai) + J^ (^0 + nai + ja^+i, ^0 + nai + (j + l)ai+i) ; 
i=o i=o 

cf. Figure |6} Integrating Dub along F^, we obtain 

1^(01 



/ Du 


= 


/ /3 + ^ 




r^|i/2 ' 


1: 


xi/2 


/,■ 



<ier/^ii/3ii2+ 



(5.14) 



using the Cauchy-Schwarz inequality. We now bound the final term in (5.14) 



Lemma 5.10. Suppose u G W ' satisfies the DOCP (5.11) and the straight cuts property 



(cf. Lemma 5.5); then 



min{|eP,#C+[a]}. 



Proof. First, we note that for any straight cut z', 



< 1- 

This follows from the fact that all 6 G -B for which z'^ = 1 can be written as 
(^ + naj,^ + naj + ctj+i) or (^ + naj, ^ + naj + aj+2) 
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Figure 6. An illustration of the definition of T^. 



for some ^ G A, ri G N and some nearest neighbour lattice direction Oj, and the definition of 
r^. Using the straight cuts property, we find that 



< 



E 



<2#C+[«]. 



5.8 



To prove the second bound, enumerate the bonds bj G F^, beginning at with the bond 
h = (^c^o + Oi)- Each bj G dC'^^ for some C^^ G C. Since u satisfies the DOCP, Corollary 
z^ 7^ implies that if z"^ connects to C^ G C^ [a] , then 

hop2(C+, C^O < hop2(Co, C'^) < hop2(Co, C'°) + 2j, 



where we have repeatedly applied the triangle inequality (2.3). Further application of the 
triangle inequality implies that the sequence of sets 

{CeC\ hop2(C, C'^) < hop2(Co, C^") + 2j} 

for j = 1, . . . , |F^| is increasing, and it is easy to see that 

#{C G C I hop2(C, C'^) < hop2(Co, C^°) + 2j}< f. 



Since by (|2.4|) each cell can contain at most one dislocation core, we must therefore have 

<|rd'<ier- 



D 



Next, we show that the number of cores can be bounded in terms of /3. This is intuitive, 
since, as discussed in [28j, each core stores a positive amount of elastic energy. Let C G C^[/3], 
then Jensen's inequality implies 

2 



|/3p> 



ac 



/3 



ac 



1 
3' 



Hence, we obtain 



#C= 



<3 



(5.15) 
(5.16) 



Proof of Lemma 5.9\ We now combine Lemma 5.10, Lemma 4.4 and (5.16) to estimate 

\{5s{o),u)\<j2\m)MO\, 



CeA 



< 



E{ii/5||2i^i"'^'+^i^(ii/^ii2'i^r)ier 



CeA 
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\i\>m\ 



lei 



-3 



We note that the |^| ^/^ G ^""^(A), so that the first term is bounded above by C||/3||2- The second 
term splits into 

CGA l«l<ll^l|2 

Straightforward radial estimates yield the bounds 

E l^r'^ll/3||2 and E 

ICI<I!/3||2 |C|>i|/3||2 

Combining the previous estimates, we obtain the stated result. D 



\^\-'< 



-1 

2 • 



5.5. Estimating ^jjgg2;b'?/''(afe). In this section, we estimate the second group in (5.3). Unlike 
in our previous estimates, which are generic, we now resort to precise quantitative bounds based 
heavily on our assumption {ip5). 



Lemma 5.11. Suppose that u G 



'1,2 



satisfies the DOCP (5.11), then for any cut z' 



connecting to a dipole C^ G C+[a] to C G C [a], we have the lower hound: 

E, , X „. , arcsinh(2/A/3) , ^ -, 

IT 



(5.17) 



beB 



where Cq> is independent of u. 



Remark 5.12. 



We have no reason to believe that the estimate (|5.17p is sharp. It is sufficient 

D 



for our purpose, due only to the fairly strong technical assumption (?/'5) made in ^3.1 



The complete proof of Lemma 5.11 is given in Appendix |X| but since this is a crucial part 
of our analysis, we provide a brief sketch. 



Sketch of the proof of Lemma 5.11 A crucial consequence of the DOCPis (5.13), since this 
inequality says that the dipoles described by each z^ satisfy 

either: hop2(Co,C^) < hop2(Co,C~); 

or: hop2(Co, C-) < hop2(Co, C+), but hop,{C+, C') < hop^{Co, C"). 

Since there is a positive dislocation core present in Co, dipoles in the first category give a 
positive contribution to the sum, since the repulsive force between C~^ and Co dominates. 

In the second case, the attractive forces between Co and C~ dominate, hence these dipoles 
give a negative contribution to the energy. Requiring the straight cuts property derived in 
Lemma 15.51 allows us to obtain estimates on the terms in the sums 

E^'(«.)^r 

fees 



using (4.9 ), and estimating sums in terms of integrals that can be evaluated explicitly. Thus, we 
obtain explicit bounds for various different cases in terms of hop2(Co,C~) and hop2(C"^, C~). 
These estimates are logarithmic, and hence turn out to essentially depend upon the ratio 

hop2(C+,C 



hop2(Co,C- 



< 1; 



the final result is the bound stated in (5.17). 



D 



We have now collected all estimates required to obtain a coercivity result. Although we state 
the result for general u G #'^'^, we will only require it later for u satisfying the DOCP. 



Theorem 5.13. Let u G 



'^''^ , then there exists a G [Dy + Du] and (5 = a — a, such that 
^(m)>ci||/3||2 + c2#C±[/3]-C3, (5.1^ 



where Cj > are independent of u. 
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Remark 5.14- We note that since we will show ||/3||2 ^ #C^[/3], we could write (5.18) more 
concisely as 



£(u) > ci 



-C3. 



D 



Proof. The bound is clearly invariant under the vertical shift and horizontal shift and rotation 



transformations we applied in ^5.2 and ^5.3 Without loss of generality, we may therefore 



assume that u satisfies the DOCP (5.11). Moreover, according to Remark 2.1 we can choose 



a E [Dy + Du] in such a way that a^ = on any bond that lies on the intersection between 
two cores. 

Summing (|5.17) over m, we find that 



J2z,^'{a,)>-r{0)^^^^^^^^^i^Cm-Co Y. hop,(C,Co)-\ (5.19) 

bes Cec-IB] 



for some Cq > 0. 



The second group in (5.19) can be estimated using the fact that each cell can contain no 



more that 1 dislocation core, so that for any S > 0, there is a constant Cs such that 

J2 hop,iC,Co)-'>-Si^C^[/3]-Cs. 



C£C- 



Bringing together (5.3), (5.9) and (5.19), we have that for arbitrary e > and 6 > 



£iu)>{ir{o)-e)m\i-(rio) 



arcsinh(2/v^) 



2tt 



+ 5 #c±H-a,^ 



we obtain from (5.15) that 
Since 



Since we assumed that «;, = on any bond that is adjacent to two cores (cf. Remark 2.1), 

2>|#C^[/3]. 



1 arcsinh(2/\/3) 

- > ^ ^ ^ ^ 0.314, 

3 71 



the result follows by taking e and 6 small enough. 



D 



5.6. Existence of Minimisers of S. With the coercivity result of Theorem |5.13| in place, we 
are now in a position to apply the Direct Method and establish existence of a minimiser of £ 
in #1-2. 

Take a sequence u" G W^'"^ such that 

£(m") ^ inf S{u). 



Referring back to ^5.3, we may assume that m" satisfies the DOCP. Let a" G [Dy + Du^' 



satisfy the condition of Theorem 5.13 and /3" := a" — d. Theorem 5.13 then implies that f3"' 
has a weakly convergent subsequence in £'^{B). In the next lemma, we also obtain convergence 
of Dm". 

Lemma 5.15. Suppose m" G W^'"^ is a minimising sequence for which each u"" satisfies 
the DOCP (5.11). Then there exists a subsequence which converges weakly in W^'"^, and the 



corresponding z^ converges weakly in any (P^B) with 1 < p <2. 



Proof. Corollary 5.13 implies that, selecting a subsequence of u" (not relabelled), we may 



assume that (5^ ^ (5 weakly in l'^{B) and that M := #C"'"[/3"] is constant along the sequence. 
Let DvP" = P^ + z^ . Lemma [573 implies that each z" can be decomposed into 

M 

771=1 
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Let B be any finite sum of positively-oriented cells. Since weak convergence in £^(i3) implies 
pointwise convergence, it follows that 



^" = - / z" ^N eZ as n ^ oo. 

dB JdB 

We enumerate the cores C^'"^ g C^[a"], n G N, m = 1, . . . , 2M + 1. Let A^bdd be the set of 
indices of cores that remain at a bounded distance from the origin, that is. 



M 



bdd := s m e {1, . . . ,2M + 1} sup dc-n,™ < +00 \. 



Since the core centres x*^"'™ with m e A^bdd can only take a finite number of positions, we can 
extract a further subsequence (not relabelled) so that they are constant. 
We therefore observe that 

A:=J2<^^[(^]= Yl *^""' for all ri. 

meA^bdd 

Then, for any finite sum of cells B, B D A, we have 

lim / /3" = lim / /3" = / (3 =: N e Z. 

"^~ JdB "^°° JdA JdA 



We aim to show that A^ = 0. For all n sufficiently large, applying Jensen's inequality implies 



in' > ^ 



/3" 

OA 



2 



N 



OA \dA\ 

Let Sq D A be a finite sum of cells that form a convex lattice polygon, and let 

Bl = BUu{C:dCnBU^%], 

for fc G N. By considering all possible corners for a convex lattice polygon it is straightforward 
to show that \dB^\ = \dB^_^\ + 6. Since lim„_s.oo /^^o /3" = N E Z, there exist Uk such that 
IdB<? Z?"'-' = N ioT i = 1, . . . ,k, and hence, 

A: „ fc 

j=0 -^^^i' i=0 ' ^' 

Since ||/3"''||2 is bounded, we obtain that A^ = 0. We can therefore conclude that, for any finite 
sum of positively oriented cells b, B D A, 

2;"" — )■ as n — )■ 00. 

dB 

Using a concentration compactness argument, we now show that we can "group" those cores 
which diverge into sums of cells with net Burgers vector zero. 

To that end, define the lattice translation operators F^ := F^^ as in (5.9), where C^ G 
C+[/3"] and /?;;:= /3" o F;;^. Note that 

11/3" II = ll/3"ll 

1 1 ' ^^ 1 1 2 I r 1 1 2' 

SO /3^ is a bounded sequence for each m, and we can select a subsequence such that /3^ ^ Pm 
for some P^ in £'^{S) and for each m = 1, . . . , M. As above, it follows that there exists a finite 
sum of positively-oriented cells A^ which contains Co and all dislocation cores of f3m, and is 
such that for any sum of positively-oriented cells B containing A^ as a subsuni, 

- f 2"oF^= /" (31^ = forn sufficiently large. (5.20) 

JdB JdB 
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Let A^ := {F^)~^A"^, then we have shown that, for all n sufficiently large, all dislocation 
cores in /3" lie within the set 

M 

S"" := U ^-- 

m=l 

We are now ready to establish that ||;z"||i is bounded. Since each z" satisfies the DMCP, 
ll^^^lloo < M. We therefore simply need to rule out the possibility that #supp{z"} — t- oo as 
n — 7- oo (i.e., that the branch cut lengths diverge). 

Fix some m G {!,..., M} and suppose that the core C^ G C+[/3"] is connected to K]^ G 
C-[/3"] via a "cut" z"'™. We claim that K^ G A""'"". 

If this were false and K]" G A^ where ii ^ m, then there must be L" G C+[/3"], L^ G A^ , 
which is connected to another core K2 G C~[/3"] outside of A^ and outside A^. Upon iterating 
this construction, we find a series of cores K^, L", K2, L2, . . . , which must eventually repeat. 

Let L" := C^. We know, by construction of the groups A^ that hop2(-Z^", -ft'j'Yi) ~^ '^^ ^^^ 
hop2(-ft'j"', i^") is bounded as n — )■ cxd for each i. This clearly contradicts the DMCP and hence 
the DOCP. 

Hence, the claim that K^ G A"'*" follows, and this immediately implies that ||2;"||i is bounded. 

Since ||2;n||i is bounded and i^ compactly embeds into £^ for any p > 1, it follows that we 
may extract a further subsequence which weakly converges in some £^(i3) with 1 < p < 2. 
This further imples that both z" and /3" converge weakly in i'^{B), and in particular that m" 
converges weakly in #'^'^, as required. D 

We now complete the proof of our main result. 

Proof of Theorem \4.5[ Invoking Lemma 5.15 suppose that the minimising sequence u^ ^ u 
in #^^'^, which is a candidate minimiser for £, and furthermore that the corresponding z"^ 
converges weakly in (.'^{B) with 1 < p <2. It remains to show that liminf £^(m"') > £{u). Recall 
from ( |4.12[ ) that 

£{u) = ^ [ip{ab + Dub) - i){ab) - 4)'{ab)Dub] + {S£{0),u). 



Since the second term is a bounded linear functional, it is weakly continuous. Using (4.9), it 
may be shown that the linear functional L defined to be 



Liz):=J2^'i 



OibjZb 



6eB 



is in (£^(i3)) for all p < 2, and is therefore also weakly continuous along the sequence. From 
(5.4) we obtain that there exists Rq, A > such that, for db > Ro, 

tlji&b + Du'i) - ^{ab) - ^'(a6)/3r > \K? > 0- 
We can therefore apply Fatou's lemma to obtain 

lim inf S^ il){ab + Du'l) - tp{ab) - ip'{a)(3J^ >y^i^{ab + Dub) - ij{ab) - ij'{ab)/3b, 

beB beB 

and thus, in combination with the weak continuity of L and 6S{0) along the minimising se- 
quence, 

inf £{v)=\immi£{u'')>£{u)> inf £{v), 



completing the proof of Theorem 4.5 D 
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6. Conclusion 

We have presented a model which has allowed us to analyse the stability of screw dislocations 
under anti-plane deformation, and we have obtained the surprising result that single screw 
dislocations exist as globally stable states, i.e. they are global energy minimisers among all 
finite energy displacements. Further, we then showed that configurations with arbitrarily many 
dislocations of arbitrary sign are locally stable, as long as the dislocation cores are suitably 
separated, but such configurations do not appear to be globally stable. 

Our work suggests two immediate directions for future study. Firstly, our analysis relies 



crucially on the technical assumption {ip5) made in ^3.1 It would be of interest to understand 



the extent to which this assumption could be weakened, but to do so would require a more 



qualitative estimate than the quantitative one made in Lemma 5.11 , and would require a deeper 
insight into the geometry of dipole interaction. Secondly, it would be interesting to understand 
which other lattice defects are globally stable states; this seems far from clear to us at present. 
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Appendix A. Analysis of Branchcuts: Proofs 



In this appendix, we detail the proofs of various geometrical lemmas from ^5.2 



A.l. Proof of Lemma 5.3, We will prove this lemma algorithmically. First, ii Zb = every- 
where, then the result is trivial. Next, we note that /3 contains an even number of dislocation 
cores, and ^C~^[/3] = #C~[/3] < +oo since u e W^'"^. We therefore enumerate C^ G C+f/?]. 

Put Cifi = Ci- Since Ci is a positive dislocation core and Zb is integer-valued, it follows 
that at least one bond bi G dCi satisfies Zb > 0. Let Ci^i be the cell such that —bi G Ci^i. 
There are now 2 possibilities: either Ci^i G C^[/3], in which case we stop, or we can find another 
bond 62 £ dCi^i such that Zb > 0. Iterating, we obtain a (possibly infinite) sequence of cells 
Cij and bonds bj. 

We now claim that no two bonds bj = bk with j 7^ fc in this sequence, and consequently 
the sequence terminates; suppose the converse for contradiction. Let j and k be indices such 
that j < k and k — j is minimal over all pairs of indices such that bj = b^. Define a polygonal 
curve P passing through the barycentres of the cells x'"^'^, . . . ,x*"i'=. P is a simple continuous 
closed curve, since the cells Cij, . . . , Ci^k-i are distinct by definition, and Ci.j = Ci^k- Hence, 
P partitions M^ into a bounded set Q (the interior of P) and an unbounded set, M^ \ Q. Define 
u = m=p1q, taking the sign according to whether P traverses dQ in an anticlockwise or clockwise 
direction respectively. 

It is now straightforward to check that Du = Du except on the bonds bi. For each of the 
bonds 6i, 

Dub^ = Dub^ -1 = I3b^+ Zb^ - 1, 



but since Zb^ > 1, this contradicts the DMCP, (5.6). It follows that the sequence bi contains no 



two identical bonds, and as Zb has compact support, the sequence must terminate at a negative 
dislocation core. 

Define z^ = ±1 if ±b G {&i, . . . ,&„}, and zl = otherwise, and then consider iterating the 
procedure described above starting at C^, but using the criterion at each step that each bond 
in the sequence must satisfy 

zb-Y,^r>o. 

m=l 
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This leads to a sequence 1-fomis, 2;*; the resulting 1-form 

TV 
i=l 

must be identically zero. If not, then the same technique as used above shows that that either 
the DMCP is violated, or else suppjz} is infinite, and hence u ^ #'^'^. 



To complete the proof of Lemma 5.3 , we simply need to show that for 2 adjacent bonds in the 
support of any z™, dbiddbi^i is a single 0-cell. This is clear, since by definition 6j, — 6,j+i G dCm,i, 
hi ^ bi+i and \dCm,i\ = 3. 



A. 2. Proof of Lemma 5.5[ We now show that we can always choose z"^ to be made up of 2 



straight cuts. Recall that Corollary 5.4 states that ||-2;™||i = hop2(C^,C~ 

Next, define the 2-cell hop operators H^ ioT i = 1, . . .6 which act on 2-cells by taking a cell C 
to the first cell 'in the direction aj', that is, the positively-oriented cell which satisfies C 7^ C, 

It is straightforward to check this is well-defined; see Figure [7j We can represent paths in the 
dual lattice by words taken from the alphabet of operators {Hi, . . . , H^}. In general however, 
the representation is non-unique - to see this, it is clear that for C2 in Figure [7| H2C2 = HiC2- 
As with the vectors a^, we define 

for any m E^L. 

We note the following properties of the operators Hi, which may be easily checked: 

(1) The orbit of the group of operators generated by {H^, . . . , H^} acting on any cell C is 
a lattice. 

(2) The operators H^ and Hj commute for any i and j. 

(3) For any cell C and any i, one of the following is true: HiC = ifj+iC or HiC = Hi^iC. 

(4) If HiC = HjC, then HiHl'^C = HjHl"'C, and if HiC ^ HjC, then HiHf^C ^ 
HjHl^C. 

(5) Hi^^Hi = Hi+^Hi is the identity map for any i. 

(6) HiC 7^ Hi±2C for any cell C and for any i. 

Setting A^ := ||2;™||i, we claim that for any pair of cells C^, it is possible to write shortest paths 
in the dual lattice as a word of the form 

H!i-,'H^ (A.l) 

for some i and some k G {1, . . . , A^}. By the definition of the hopping operators, it is clear that 
such words represent a sequence of cells lying on the lines x*^™ + tai and x'"™ — taj+i. 

We now prove the claim: first, we show that any shortest path must be able to be written as 
a word made up of only 2 of the operators Hi and ifj+i. Suppose this is false, for contradiction. 
(6) implies that we may assume that it contains a segment which may be written as 

Hi+2HpHi_2 or Hi_2HJ"Hi^2 

for some i and some m G N. Since both cases are similar, we consider only the first. If C is the 
first cell in this subsequence, Hi_2C 7^ Hi_^C, or else HiHi^2C = HiHi^^C = C by (5), and 
there exists a shorter path. (3) therefore implies 

Hi_2C = H.^iC (A.2) 

Invoking property (3) again, either 

HiHi^2C = HiHi^iC = H^__iC or HiHi^2C = Hij^iHi^2C = C. (A. 3) 

Once more, the second case results in a contradiction; hence repeatedly using (2), 

-ni+2-ni ni_2^ — ^i rti+2^i^i-2'^ — -"j J^i+2^i-i'-^ — -"j ^i-l'^, 
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Figure 7. This illustration shows the operation of the operator Ha^ on two 
representative cells, Ci and C2. 



implying a contradiction in the case where m is odd. In the case that m is even, repeatedly 



using (2), (A. 2) and (5), we have 



r2fci 



T2k, 



TT Tj2k TT r^ u"^^ u u n u'^*^ u u n u'^^^r^ 

-ni+2-nj J^i~2<^ — -Hj -nj+2-ni-2<-^ — -Hj -ni+2-ni-l<-^ — -"j <--, 

obtaining another contradiction; it follows that it is impossible that every shortest path can be 
written as a word containing no more than 2 of the operators Hi. 

We next show that if the two operators are Hi and -ffj+2 for some i, we may rewrite the word 



in terms of Hi and Hi_i, or Hi^i and Hi. In (A. 3), we showed that for any cell C contained in 



a shortest path HiHi_2C = Hf_^C. Suppose that a shortest path is represented as a product 
of Hi^2 and H^. Then each pair HiHi^2 niay be replaced by H]i_^, and using (2) to permute 
each of these pairs to the right, we eventually obtain one of 



HTHt 



2n 



or 



Ijm TT TT2n 



In the first case, the proof is complete. In the second case, it must be that HiH,^ 
else m = 0, since otherwise (3) implies 



2n 
i-1 



H 



2n+l 



or 



TJTn TT TT2n 

-Hi_2-ni-ni-l 



rrm tt Zj2n 

-"j-2-ni+i-^j-i 



Tjm—l TT2n 
-Hi_2 -",_!, 



which is a contradiction. Hence we have proved the claim, and in fact since we have obtained 



a shortest path in the form (A.l), the lemma is proven for this particular case. 



Finally, we consider a general word made up of only the operators Hi and -ffj+i. Now consider 
a word of the form Hi^iHJ^Hi^i. If m is even, then we can generate a new word corresponding 

If m is odd, then we can write a new shortest path as 

U U U Tjm—l 
-ni+l-nj-ni+l-ni 



to a shortest path Hf^^^H"^ 



But then using (3), this must be able to be written either as 



Hi+iH. 



m—l 



or Hi+2Hi-iHi+2H, 



m—l 



where the second case results in a contradiction. It is now possible to check that this implies 
the full conclusion, since by these arguments we can always transform a general word composed 
of Hi and ifj+i in one of the form ( A.l ). 



A. 3. Proof of Lemma 5.11 , To prove this lemma, we use Lemma 5.5 to assert that z should 
be made up of straight cuts. This will allow us to make estimates for each straight segment, 
which depend upon on the orientation of each segment. We divide the lattice into 'sextants' 
by defining 

iS* := |x G M^ I 3 A > 0, /i > such that x = Xai + /iOj+i j. 
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and rings TZ^ by defining TZ^^ = 0, 7^° := C'o, and then 

n' = cios \J{c ec\ [TV-^ u n'-^) \ cios(C) n dosijv-^) ^ 0}. 

We will say a straight cut z is: 

(1) tangential if suppj^;} C 5* and the cut direction is aj+2, or equivalently suppj^;} C 
5* n TV for some r, and 

(2) radial if suppj^;} C 5* and it has direction aj or Cj+i. 

It may be checked that for any straight cut, there exist C^, C^ G C such that 

z = ±1. 

We separate the full result into 2 further lemmas, which give precise estimates for each of 
these classes of straight cuts, before combining them to complete the proof for the general case. 
By Taylor expanding il)'{ah) around we write 

Y,i'\&b)zi, = Y,'i'"^^)^^^b+\i'^^\sb){a^fz,. (A.4) 

feez? b&B 

We proceed to estimate the first terms in the summand, by estimating on the radial and 
tangential straight segments of z"^ separately. 

Lemma A.l. For a tangential cut z^'^^ with ||-2*'^'^||i = hop2(C^, C~) = I on ring TV , we have 
the following estimate: 

Proof. First, we appeal to symmetry. If 6 = (.^, .^ + flj) G iSj, then applying the refiection 

R = -{ai + aj+i) (g) (flj + aj+i) - 0^+2 ® 0^+2, 
it is straightforward to check that if b' = {RC,, RQj then 

This means that if a tangential cut crosses the line of symmetry {t(aj + flj+i) \t G M}, then 
some of the bond contributions cancel. It follows that we need only consider the case where 
all bonds in the support on 2;*^'^ lie on one side of this line of symmetry, since this is the worst 
case. For such cuts, it may be checked that / < r. 

Identifying z^^"^ with a dual lattice path, we may enumerate bk G suppj^;*^"} 'along the path', 
and elementary geometry now shows that a^^, has alternating sign as k increases. Letting 
^^ := |(r — |)(aj + Oj+i) if i is odd, ^^ := ^(r — ^){ai + flj+i) if i is even, it is straightforward 
to show that each bond b^ can be represented as one of 

{C + sai+2, r + sai+2 + tti) or {^ + saj+2, C + sai+2 + Oi+i) 
where s G |0, |, . . . , ^-y"}- Elementary trigonometry now allows us to calculate that 



27rdfe 4^" = ± 



or ± 



s + l/2\ / s 

arctan | ,^ , , , — arctan 7-— - 

, s-l/2\ / s 

arctan ,^ , _,, — arctan 7— - 
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respectively for the cases above. We therefore have that 

\l/2\ 

27r Y^ zf'ab = J2 arctan (^^^f^) - 2 arctan hpf) + arctan (^^f^) + 0(r-^), 

b t=l 

[1/2] 

or = 

t=i 

where Sq E {O, ^, . . . , ^^}, and the 0{r~^) term arises from the contribution of at most 2 



1^ 



J2 arctan (^^^^^^ J - 2 arctan (^^ + arctan (^±^|^) + 0{r-^), 



bonds we have neglected, whose contribution we estimate using (4.9). In the second case, the 
fact that — arctan is convex for positive arguments implies that the sum is bounded below by 
0. In the first case, for all t in the range of summation, 

t + so + 1/2 t + So 1 , • • 1, t + So t + So 
1^ ^11 > 1^ I and trivially t— — tt < .^ , • 

Since arctan is increasing, we obtain the lower bound 

[1/2] 






t=i 

> — arctan 



using the fact that arctan is positive and increasing for positive arguments. Finally, note that 
in the case where I > r, i.e. the tangential cut crosses the line of symmetry, and we obtain the 
same estimate but with 2r — / in place of / in the formula above, so using the definition of ^'^ 
gives the result. D 

Lemma A. 2. For a radial cut z^''^'^ such that ||z'"'^'^||i = hop2(C^, C^) = I, either \x^ \ < \x^ \ 
and 

beB 
or else \x^ \ > |x*" |, and if C^ E W , then 

Yzl^^a, > --arcsinhf— 5^/^^ ^^ -0{r-'). 



beB 



Proof. First, we enumerate the bonds in bk E {b E B \ zl^'^ > 0}, beginning with the bond for 
which db is smallest, and proceeding outwards along the cut. Elementary geometry shows that 
the terms zl'^'^atk are all positive in the case where bi is in one of the directions Oj+i, aj+2 or 
Oj+a, which corresponds to having \x'" \ < \x'" |; this immediately provides the first bound. 



In the second case (4.9) implies 

^b, > 



2'Kdhf, 

Without loss of generality, we assume the cut direction is a^, the case with direction aj+i being 

similar. There are now two cases: bi is either in the direction aj_i, or aj_2. Further elementary 

geometry allows us to conclude that in the first case, d^^ = \x\ for some x E S\ and in the 

second, di,^ > \x\ with x = x^ E clos(5*). In either case, di,^ satisfies the same lower bound as 

dbj^, and further, we have that 

db2n^i^db2n ^ \x + nai\. 

Noting that as x G clos(5*), it follows that a, ■ x > and 

1 1 1 

< 



\x + tai\ ^\x\'^ + 2tai-x + t^ ^\x\^ + t^ 



SCREW DISLOCATION UNDER ANTI-PLANE DEFORMATION 



31 



which is a decreasing function of t, so we estimate 
I . R/21 



k=l i=0 



1 



> 



\X 



nai 



1 

\x\ 



\l/2\ 



V~\ 



ds 



xr + s^ 



Evaluating the integral, and noting further that \x\> ^(r — 2/3), we obtain the conclusion. D 



Vli 



We now combine the estimates of Lemma lA.ll and Lemma IA.2I to obtain an estimate for a 
general cut z™ made up of two straight segments. As we showed in Lemma |5.5[ each z™ is 



made up of at most 2 straight segments. It may be checked that each of these segments is 
either purely radial, purely tangential, or changes from tangential to radial part way along its 
length, with one bond which crosses dS^ . All possible cuts satisfying the DMCP and made 
up of 2 straight segments can therefore be decomposed as either 

(1) a tangential cut and 2 radial cuts or 

(2) a tangential cut, a radial cut and tangential cut, 

where any of these segments could possibly have length 0, and neglecting the extra bonds 
mentioned above for now. Recalling the result of Corollary 5.8[ 

hop2(Co,CJ>hop2(Co,C+). 

Consider the first case, letting the two radial segments be of length li and I2 respectively, and 

*^'" I and the radial cuts are of non-zero 






< \x 



the tangential segment of length I — li — I2. li \x 

length, then the radial segments have the trivial lower bound, by Lemma A. 2 In the worst 

case, where I — h — I2 = r, we have the bound 



Otherwise, 



E 

r+ I 



z?a,> 



27r 



arctan 



r — 



2/3)v^; 



0{r- 



1 
12 



0{r- 



\x 



> \x' 



Urr 



SO applying Corollary 



5.8 



we have that for C E TV 
2r > hop2(Co, C") > hop2(C+, C) =h + h. 
Therefore, applying Lemma [A. 1 and Lemma A. 2 



(A.5) 



beB 



27r 



arctan 



2[(/-/i-/2)/2j +1 



(r - 2/3) v^ 



+ 2 arcsinh 



2 arcsinh ( 



( mm 

\V3{r-2/3) 

2\l2/2] 



VV3(r+[/i/2j -2/3) 



0(r- 



By ignoring the floor functions, it is possible to check that under the bound (A.5), the function 



in parentheses is increasing in both li and I2 if r is suitably large; we therefore have that the 
maximum must occur when I = li + I2. Hence putting I2 = I — h, we have 



> Zbab> arcsmh —= 



baB 



2/3) 



+ arcsinh ( 



2\{l-h)/2\ 



Vy3(r+[/i/2j -2/3) 



Oir- 



Dropping the floor and ceiling functions, this estimate is convex in /i. The worst cases are 
therefore /i = or /i = / and 2[//2] = 2r + 1, giving the value 



1 / 2r + 1 \ 

> Zhah> arcsinh! —pz | + Oir^^] 

ti -^ \V3{r-2/3)) 



arcsinh (2/\/3) 



+ 0(r 



-l^ 



TT 



In the case where we have tangential, radial and tangential segments, similar arguments show 
that, once more, the worst possible bound arises in the case where the cut is purely radial, 
giving the same lower bound. 



SCREW DISLOCATION UNDER ANTI-PLANE DEFORMATION 32 



We now return to the lower order terms in (A.4). Then by crudely estimating 



5^ i^(^) (.,)(«,)%- <j:d^'<i 



beB 



r 

dh>r 



and noting that the worst case bounds always have occur when hop2(Co, C^) ~ r, we have 

E.,. „, , arcsinhf2/-\/3) 
^^'(Afc) > -^"(0) ^^^^ - cohop2(Co,C; 



[1 

[2; 

[3 

[4 

[s; 

[6 

[7; 

[9 

[lo; 
[11 
[12; 

[13 
[14 

[is: 

[16 

[1?: 
[is: 

[19 
[20 

[21 
[22 

[23 
[24 

[25 

[26 
[27 



TT 



-1 

'ml 
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